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TECHNICAL MEMORANDUM 


DOPPLER-CANCELLED RESPONSE TO VLF 
GRAVITATIONAL WAVES 

I. INTRODUCTION 


On several occasions [1,2] it has been pointed out that Doppler 
tracking of a distant spacecraft can , under certain conditions , provide 
a valid technique in the search for Very Low Frequency (VLF) gravita- 
tional radiation. The basic system consists of the Earth and a distant 
spacecraft considered as free masse.; communicating through microwaves: 
a monochromatic electromagnetic v» ve is emitted at a station, relayed by 
the spacecraft, and, finally, received at a station where its frequency is 
compared to that of a reference oscillator. The resulting frequency shift 
is primarily caused by the Doppler effect. Gravitational waves have been 
suggested among the potential causes of Doppler residuals [ 3] . More 
precisely, it has been pointed out [4] that those plane gravitational waves, 
for which the product of their frequency by the round-trip light time is 
a large integer, would produce preferred peaks in the spectra of the 
Doppler residuals. This result is based on a "two-way" configuration in 
which the light signal interacting with the gravitational wave follows a 
round trip [5). Similar spectral properties can, in general, be expected 
for combinations of microwave links other than the two-way type just 
mentioned, but the response to gravitational waves will change 
accordingly. 

Recently, NASA and the Smithsonian Astrophysical Observatory 
(SAO) [6] have accomplished a verification of an important quantitative 
prediction of the principle of equivalence. The experiment compared the 
frequency of a space-borne atomic clock with that of an identical clock 
at an Earth station by using one-way ( spacecraft - Earth ) and two-way 
(Earth- spacecraft -Earth) Doppler information. This experiment in par- 
ticular demonstrated that under certain conditions a combination of one- 
and two-way microwave links provides, in real time, a beat signal which 
is independent of the otherwise dominating Doppler effect but does depend 
on other weaker effects. 

The system which removes the Doppler effect has been called the 
Doppler Cancelling System, and the Doppler-cancelled signal is the 
corresponding observable [ 7] . 

Calling ^2w accumu ^ ate( ^ cycles in the one-way down 

link and two-way link, respectively, the phase of the Doppler- cancelled 
beat signal is defined by 



( 1 ) 


s * * s *lw • 5 s *»w 


and its fractional frequency by 




( 2 ) 


where f is the carrier frequency and Af/f) lw and Af/f) 2w are the fre- 
quency ahifta in the one-way and two-way signala. 

By definition, the Doppler cancelling technique will cancel every 
frequency shift which ia doubled when the path length of the radio wave 
is doubled. For plane, long periodic gravitational waves, which may be 
considered as inducing time-dependent perturbations in the time of prop- 
agation of the radio signal, the frequency shift caused by the space-time 
curvature during the round trip is not twice the frequency shift corres- 
ponding to the one-way path, in general. 1 Thus, the Doppler cancelling 
technique is expected to respond to gravitational waves of sufficiently 
low frequency for geometric optics to apply but high enough to resonate 
the Earth- spacecraft system and avoid local measurements. 

The purpose of this work is to identify the response of thi Doppler- 
cancelled observable to pulses and continuous gravitational waves. Section 
II treats the effect of the curvature of a plane gravitational wave space- 
time on signal propagation. The induced frequency shift on the light 
signal is computed by means of a new approach which uses the well- 
known [8,9] analogy between a metric gravitational field and a medium 
with an equivalent index of refraction. The result is compared with that 
of a previous calculation [4] based on a different approach, and they are 
found to be identical. In Section III the response of the Doppler Cancell- 
ing System to bursts of gravitational waves and continuous radiation is 
identified. Response functions and spectra are constructed by vising the 
mathematical methods recently proposed [10] to deal with gravitational 
wave effects on light propagation. 


II. EFFECT OF A GRAVITATIONAL WAVE 
ON SIGNAL PROPAGATION 


A plane gravitational wave is described, in the "transverse, trace- 
less" gauge [11], by a space-time with a Lorentz metric 

T. This holds provided that the wavelength of the gravitational wave is 
not too large in comparison to the Earth- spacecraft distance. Other- 
wise, the locality of the measurement implies zero response for every 
Doppler technique. 


2 


( 3 ) 


ds 2 * n w <** a d* b - (dx 2 - dy 2 ) + 2A„ dx dy , 
ao + x 


where the approximately inertial frame of reference has been oriented so 
that the gravitational wave travels along the z-axis. The amplitudes A + , 

A x of the two independent states of polarization are a retarded solution 

of the wave equation in Minkowski space-time 


Aj = Aj (ct - z) , i * +, x 


(4) 


Pointwise they satisfy the weak field condition 



(5) 


Perturbations of the Minkowski metric in general cause time delays in the 
propagation of light signals. For a wave-like perturbation , as in equation 
(3), the variations in accumulated phase and the frequency shift induced 
on a monochromatic electromagnetic wave can be computed using two 
approaches. One technique [4] consists of computing the frequency shift 
by integrating the equation of parallel transport of tetrad components of 
a null vector, the four-dimensional wave vector, defined by the metric 
(3). A different technique is used here. It consists of treating rela- 
tivistic effects on a null path by geometrical optics in a three-dimensional 
space with an equivalent index of refraction. This procedure gives 
directly the time delay due to propagation and the corresponding accumu- 
lated cycles. Time differentiation gives the number of cycles counted in 
the unit interval of time, and the frequency shift is obtained. 

Work to linear order in A and 1/c. Consider that a light signal is 
exchanged by an Earth station and a distant spacecraft and define the 
cosine directions of the geocentric position vector of the spacecraft 


u = % , s = Hf • Y * a£ : dr • <o * 2 ♦ “y 2 + d*V . 

( 6 ) 


Note that (Fig. 1) 


8 2 - y 2 = (1 - u 2 ) cos 24> ; 20y = (1 * uh sin 2<J> , (7) 


<J> being the longitude of the spacecraft. 
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t 



Figure 1. Gravitational wave and Earth-spacecraft geometry. 
The path of a light signal in space time satisfies 


d» 2 » 0 . (8) 


From equation (3) the time delay between emission and reception of the 
light signal is formally given by 


c(t r - t e ) = / r |^1 + 1 ' £ (A + cos 2$ - A x sin 2$)J dr . (9) 


The integrand in equation (9) may be considered as an equivalent index 
of refraction. The integral must be extended from the Earth station, 
with coordinates (r^, tj), to the spacecraft, with coordinates (r 2 . tj) , 

and back to the receiving station, with coordinates (r 3 , tj). We adapt 

the origin of coordinates so that r ^ = r 3 = 0. Then, the parametric 

equations of the forward path are, in first approximation, 


4 


X = 0P 


t 


y 8 yr 


t s pr 


( 10 ) 


°I x i x 2 » 0 < y < y 2 . 0£*£*2 

which imply 

Aj(ct - z) = Aj[r(l - m)J (11) 

along the forward path. For simplicity , we set 

l * r 2 ; h = cos 2$A + - sin 2$ A x . (12) 

Evaluation of equation (9) along the approximately null trajectory (10) 
gives 

l 

c(t 2 - t x ) = i + i-J-H / ii(u) du , (13) 


where u=r(l-u)isan outgoing null coordinate. The parametric 
equations of the backwards path are, in the same approximation, 

x - x 2 = -0(r * A) , y - y 2 = -y(r - 1) , z - z 2 * -y(r 
which imply 

Aj (ct - z) * Aj (r (1 + y)l , i = +, x 

along the backwards path. Thus, 

l 

e(t 3 - t 2 ) = 4 - J ii(v) dv . 


- I ) 
(14) 


(15) 


(16) 
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where v = r(l + y) is an ingoing, null coordinate. From equationa (13) 
and (16) it followa that the total time delay for a round trip is, to linear 
order in the A's and 1/c, 


c(t 3 - t x ) ■ 24 ♦ J |^(1 + u) 

We define the emitted , transmitted , 
of cycles dN per unit of time, 2 


l l 

f h(u> du + (1 ~ u) / h(v) dv 
9 0 

(17) 

and received frequency as the number 



dN , .dN , dN 

i a r x • f 2 -ar 2 * f 3 a ar 3 


(18) 


respectively. Equations (13) and (16) give the frequency shift in the 
upward and downward paths, respectively. 


Af\ . dt l _ . _ . I . 1 ♦ u ' . y . 

tJLp *»r ‘ ‘V 


dt. 


Af\ . “ v 2 


'down ^3 


(h < 


l 2 ‘ * 


where h^= h (ctj - Zj). The Doppler-cancelled frequency shift is 

A f _ 1 / A f I Af | \ 1 + V 1 , 1 v 

T ‘ H T Lm 'T lup) ' ’ h > + * 


(19) 


M h. + i h„ - h, , (20) 


while the two-way frequency shift corresponding to the round trip is 


1 /2w /up /down 


» l+j a; i - u : 21 

5 h l uh 2 ~T~ h 3 * c 


( 21 ) 


T. These definitions do not account for the time scales of the Earth 
' and on-board clocks being offset by a nonzero relative gravitational 
field (redshift). This timing effect must be considered separately. 
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Equation (19) coincide* with the one derived by integrating the equation 
of parallel transport of the electromagnetic wave vector, i.e. , the null 
geodesic equation (4). The present approach does not make explicit use 
of the geodesic equation . This equation is, however, automatically 
satisfied by virtue of equation (8) and the fact that, in geometric optics, 
light rays are normal to surfaces of constant phase. 


III. RESPONSE TO BURSTS AND CONTINUOUS 
GRAVITATIONAL WAVES 


From equation (20) it follows that the response of the Doppler- 
cancelled observable to a sufficiently short (< t/c) burst of gravitational 
waves consists - as for equation (21) - of a threefold repetition of the 
waveform, though with somewhat different weighting coefficients. In 
particular, for orthogonally incident gravitational waves (u 8 0), the 
three-pulse structure is preserved. 

Helling* [ 10] has painted out that an optimum technique for the 
detection of a signature of known form in data corrupted by white noise 
is the matched filter. He has further shown how to construct one 
appropriate for the two-way Doppler observable. The same method can 
be used for the Doppler-cancelled observable as follows: A matched 
filter [12] is operationally defined by its transfer function 


-iu>T 

F(u>) « a S%) e m 


( 22 ) 


where a is a constant which can be precisely defined in terms of the 
noise power, S*(u>) is the complex conjugate of the Fourier transform of 
the theoretical signature, and t is the time of maximal signal-to-noise 
ratio. We define 11 


Tj 8 21 /c , t 2 8 1(1 + u)/c , Tj 8 0 


(23) 


1 + li 1 1 - u 

*1 8 - -J-K , a 2 8 2 . a 3 8 - 


(24) 


and 


y( t) » 2 “j - 
i*i 


V 


(25) 


Equation (25) shows that if a fsatura is identified in tha Doppler raskluals 
at the tins x, similar features should ba recognizable at previous instants 
Xj and ig. We assume that y(x) and 6(x) are of integrable square over 

tha real Una and denote by y(w), h(u>) tha respective Fourier transforms. 
Thus, tha gravitational wave signature is, in the frequency domain, 


3 

y(w) » h(w) J 

J*1 



(26) 


The corresponding matched filter is, using equation (22), 


F(u>) a a h*(u) e ^ 2 Sj e iwT * (27) 

Jal 

in the frequency domain and 

00 

F(x) * a J h(x') g (x - x') dx' (28) 


in the time domain, where 
3 

g(T) a J 5 (t - Ct^i, - 1j>) . (29) 

ial 

We put = Tj. From equations (28) and (29), it follows that the filter- 
ing procedure consists in matching the input signal to the theoretical 
signature run backwards from t * x m to t = Xj. 

According to some astrophysical models [13], pulses coming from 
different events may combine stochastically, producing a signal with a 
higher amplitude than in the case of a pulse from a single event. It is, 
therefore, useful to average the filter - or, equivalently, g(x) - ever 
the angle of incidence. The result ia 
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(30) 


1 

W T) ‘I / « <T) 

-1 


* - \ 1 6 (t) «■ 6 (t - 22/0) ♦ 


c/42 0 < t < 22/c 

0 elsewhere 


The corresponding (liter is the convolution of equation (30) with the 
gravitational waveform h(t). The second pulse repetition, whose position 
in time depends on v>, has been averaged out. The same happens for 
the two* way observable, which has, in fact, the same time sequence 
(equation (23)] as the Doppler-cancelled observable. What distinguishes 
equation (30) from the corresponding filter of the two-way system 
(equation (37) in Reference 10) 


t ' 1 a C 0 < t < 22/c 

<2/c) z 

» 

0 elsewhere 

(31) 

is the sign of the correlation between the spikes at t ■ 0 and t = 2l/c: 
this is positive in equation (30), negative in equation (31). 

The Doppler-cancelled response to continuous gravitational waves 
is now considered. The autocorrelation function 


*ave.2w (t) ■««>-««- 22/0 ♦ 


R(t) 


/ y(r*) y(T -t’ ) dt' = ^ / 


|y(u>) | 2 e iwT dw (32) 


is, for the process (equation (25)), 

• 3 

R(x) * ^ / |h(d>) | 2 Y ®j ®* c 0 ® • V e iwT dm 

-» j .2=1 

Ob 

* Hr? / !**(“) I 2 (1 + COS 2 m2 * u 2 sin 2 mi 


- 2 cos m2 cos m2u 

- 2u sin m2 sin mlu) e iu>T dm . (33) 
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The corresponding spectral density (one-sided) 


S (u) ) = |y(u) | 


2 


0 < (I) < 00 


( 34 ) 


is 

S(oj) = | |h(u))| 2 (1 + cos 2 itiS, + y 2 sin 2 u>A - 2 cos ujA cos wAy 

- 2y sin <*>A sin u>Ay) . (35) 


For gravitational waves orthogonal to the line of sight, 


S(b>) = |h(w) | 2 sin 4 (u = 0) (36) 


This spectrum has exactly twice the periodicity of the corresponding 
two-way Doppler spectrum [equation (24) of Reference 4] 


S 2w (u) = l h(u,) i 2 sin2 = 0) 


(37) 


Averaging equation (35) over y gives 


S ave <“> 


n ' h< “> 


( 


1 + i COS 2(dA 


5(uA) 


. 3 cos 2uiV 

I 5 , 0x2 . 
(u»A) ) 


(38) 


Equation (38) is the spectrum of the Doppler -cancelled response function 
to an isotropic background of random gravitational waves. 

In Figure 2 the autocorrelations of the Doppler- cancelled response 
(left) and of the two-way Doppler (right) are compared in the case of a 
Hat spectrum [h(co) = const.] at angles of incidence 6 = 30°, 60°, 90°, 
and isotropic. The analytic expression of the two-way spectrum is [10] 


2 2 2 2 2 

Sgw^) = |h(<*>) | (y + sin u>A + y cos u A - 2y sin ouA sin wAy 


2y cos u)A cos u)Ay ) 


(39) 
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T 


Figure 2. Three-way Doppler -cancelled (left) and two-way Doppler 
(right) autocorrelation function (one-sided). A constant spectrum 
is assumed in both cases. The angles of incidence are (a) 30°, 
(b) 60°, (c) 90°, (d) isotropic. The numerical values on 
> the y-axes are defined by the number of points (128) 

used in the calculation. 
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3 


-1 


for an angle of Incidence 6 = cos u and 


■ l h< “> |2 (* • 3 “* *“* - ^7 


cos 2uJ l . 2 sin 2wA 

t + . 


(wA) 4 


(<*>A) 4 


) 


ram 


for averaged angle of incidence. In Figure 2 it is interesting to note 
that at 0 = 90° the Doppler-cancelled autocorrelation has a peak at t * 
lie. All the intermediate features disappear by averaging over 0. The 
peak at t = 2A/c is 1/10, in the Doppler- cancelled case, and -1/6, in the 
two-way case, of that at t = 0, as is evident from the form of equations 
(38) and (40) in the limit of large values for ail. 

The numerical calculation evaluated equations (35) and (39) at 
selected values of y , and equations (38) and (40), each at 128 equally 
spaced points, with eight points in each of the intervals 2nir < id A < 
2(n+l)w, n = 0, 1, 2, ... 15. A Fast Fourier Transform routine was 
then used to obtain the corresponding autocorrelation function [equation 
(32)]. 


IV. CONCLUSION 


As has been pointed out [14,15], the use of multiple time- 
correlated radio signals offers interesting possibilities not only for redshift 
experiments but also in the search for VLF gravitational waves. In par- 
ticular, the same Doppler Cancelling System developed for the NASA-SAO 
redshift experiment responds to VLF gravitational waves in a precise 
manner. The form of the response function is identical to that of the 
two-way Doppler configuration as far as the timing is concerned but is 
otherwise different. For example, there is always a positive correlation 
in the Doppler-cancelled response function at two events separated by a 
round-trip light-time. Further, a three-peak feature is maintained in the 
case of gravitational waves propagating at a right angle with the line of 
sight to the spacecraft. This is a most interesting case, because at 
this angle of incidence the transverse nature of general relativistic 
gravitational waves would maximize the strength of the coupling with the 
Earth- spacecraft system. 

An analysis of the influence , in the Doppler Cancelling System , of 
other timing and propagation effects may indicate whether the features 
of the Doppler-cancelled response function to VLF gravitational waves 
could be decoupled from those of at least some of the remaining noise 
sources. 
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